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■ Two mechanisms are examined for hard double pomeron exchange dijet 

production, the factorized model of Ingelman-Schlein and the nonfactorized 

cn \ model of lossless jet production which exhibits the Collins-Frankfurt-Strikman 

mechanism. Comparison between these two mechanisms are made of the total 

CNj ' cross section, E'y-spectra, and mean rapidity spectra. For both mechanisms, 

several specific models are examined with the cuts of CDF, D0 and repre- 
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^^ . sentative cuts of LHC. Distinct qualitative differences are predicted by the 

^^ , two mechanisms for the CDF y+-spectra and for the -ET-spectra for all three 

experimental cuts. The preliminary CDF and D0 experimental data for this 
On ' process are interpreted in terms of these two mechanisms. The y+-spectra of 

the CDF data is suggestive of domination by the factorized Ingelman-Schlein 






mechanism, whereas the D0 data shows no greater preference for either mech- 
anism. An inconsistency is found amongst all the theoretical models in at- 



Q , tempting to explain the ratio of the cross sections given by the data from 

1^ ' these two experiments. 
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I. INTRODUCTION 

In diffractive hard scattering, the incident hadron in e-p collisions and one or both 
hadrons in pp collisions participate in a hard interaction involving very large momentum 
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transfer, but nevertheless the respective hadrons emerge with small transverse momenta 
and a loss of small fractions of their longitudinal momenta. For such diffractive hard pro- 
cesses, first comes a question of pure semantics, whether or not to say the diffractive proton 
"exchanged a pomeron". Only one Pomeron has entitled historical rights to this name, and 
that is the Pomeron of soft Regge physics ||1]J^ (also sometimes called the soft Pomeron). 
Reference to pomeron in any other case exploits this established trademark as a mnemonic 
for describing some portion of the process in which a strong interaction scattering occurred 
that involved the exchange of no quantum numbers except angular momentum. In our dis- 
cussion of diffractive hard scattering, we will use the lower case pomeron in reference to a 
process in which one or both incoming hadrons diffracts into the final state along with a 
hard process. On the other hand, the upper case Pomeron will be reserved for the vacuum 
exchange trajectory of soft Regge physics ^^■ 

There is general belief that properties of the Pomeron reflect in the pomeron of diffractive 
hard scattering, although it is a central research question to identify the specifics. Spacetime 
arguments generically suggest that hard events are well localized in space and time. Thus it 
is expected that in a diffractive hard process, the diffractive hadrons undergo effects similar 
to what they would encounter in a high-energy elastic scattering. As such, diffractive hard 
physics is expected to involve long-time, long distance, thus nonperturbative, physics. Nev- 
ertheless, that hard processes can occur intermittent to the diffractive scattering indicates 
that diffractive hadronic physics, via the pomeron, also possesses perturbative properties 
that can be explained through perturbative QCD. 

A primary goal of diffractive hard scattering physics is to unify the QCD picture of the 
pomeron with the phenomenological Regge physics description (for a review of Regge phe- 
nomenology applied to diffractive physics please see 0,H). Hard double pomeron exchange 
(DPE) processes are useful in addressing this question, since it turns out the QCD and Regge 
physics description of these processes have some distinct qualitative differences, which are 
best expressed in the context of hard factorization. 

Recall, for a hard scattering factorized process, the effect of the two incoming particles 
act independently on the hard event |^J^ . The basic Regge physics motivated model of hard 
diffractive processes is the Ingelman-Schlein modelQ 0, and this model assumes hard fac- 
torization. In their model, diffractive scattering is attributed to the exchange of a pomeron, 
which operationally is defined as a colorless object with vacuum quantum numbers. Their 
model treats the pomeron like a real particle and so considers, for example, that a diffrac- 
tive electron-proton collision is due to an electron-pomeron collision and that a diffractive 
hadron-hadron collision is due to a proton-pomeron collision for single-sided diffraction and 
pomeron-pomeron collision for double diffraction. 

For diffractive deep inelastic scattering, basic ideas of hard factorization were outlined 



^Their model was motivated by a prior and seminal diffractive hard scattering experiment by the 
UA4 S and subsequently the ideas of their model were first studied by a UA8 experiment [^. 



Some other theoretical works at around the same time as this model also had similar ideas |1C,11| 



and diffractive parton distribution functions were defined in [0,|T^]0. A proof of factoriza- 
tion for diffractive DIS was given in |15|. For hard diffraction in pure hadronic collisions, 
Collins, Frankfurt and Strikman (CFS) |T6| have demonstrated a counterargument to hard 
factorization. The CFS mechanism is a leading twist effect in which all the momentum 
lost by the diffractive hadrons goes into the hard event. An important feature about the 
CFS mechanism is that it requires the color flow properties of QCD in an essential way. In 
general, the presence of color in QCD implies pomeron exchange in simplest form is a two 
gluon exchange process [p!7| , p!8| . Necessarily, the simplest model of the pomeron must involve 



at least two partons in order to be color singlet. The two gluon pomeron model has a key 
property for any pure hadron initiated reaction, which is a realization of the CFS mecha- 
nism. Consider the hard DPE process /11/12 — > /11/12 + hard, where /ii,/i2 are the colliding 
hadrons. The two gluons exchanged by hi are not both obliged to enter the hard event. 
Instead, one gluon may attach to /i2- In this case, the two incoming hadrons no longer act 
independently in inducing the hard event. By definition of hard factorization [§],0, such a 
process is nonfactorizing. 

This mechanism was identified earlier by Frankfurt and Strikman |T^. They originally 
referred to the nonfactorized pomeron of CFS as the coherent pomeron. Subsequently, the 
UA8 presented results [^ in which up to 30 percent of the dijet events in single-sided 



diffraction could be associated with the coherent pomeron, which they in turn named the 
superhard pomeron. With the hindsight of the UA8 experiment and the ideas of CFS, in 
H^ the CFS mechanism was applied to a toy quantum field theory model of diffractive dijet 



photoproduction, in which the pomeron was represented by two gluon exchange. This work 
in turn, in turn, named the nonfactorizing, alias superhard, alias coherent pomeron process 
as lossless diffractive hard scattering to emphasize the efficient transfer of the pomeron 
momentum to the hard process. 

The CFS mechanism has been developed for hard DPE in pp collisions^ for quark jets in 
2T| and gluon jets in [^]. The gluon jet process was shown in ||2^ to dominate the quark 



jet process by several orders of magnitude. 

The purpose of this paper is to examine for the DPE dijet process, general differences 
between the factorized pomeron model of Ingelman and Schlein, F(IS)DPE ||^, and the non- 
factorized pomeron model of lossless jet production of Berera and Collins p2|, N(L)DPE. 
Our notation specifies in the context of hard factorization whether the process is factorizable, 
F, or nonfactorizable, N, and in parenthesis gives the particular type of process. The latter 
specification is necessary since there are several different types of factorizable and nonfactor- 
izable processes. Detailed discussions about this point are in [Il6| , p!2| , pr3| , p2| , p!5| , p5| , p6| . As one 



example, factorized processes first have a basic distinction between simple hard factorization 



'Closely related to diffractive parton distribution functions are fracture functions |14] 
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The first nonfactorizing DPE two gluon model was developed before CFS for Higgs [23| and 



heavy quark |23] production. Although the nonfactorizing mechanism is the same as that of CFS 



[iql and |12,21,22|, these earlier papers did not recognize the full consequences of nonfactorization 



to the extent done by CFS. 



and the more specific Regge factorization [0. In particular, tlie factorized Ingelman-Sclilein 
DPE model also is Regge factorized^. 

For nonfactorization, one example outside of the CFS mechanism is the "flux renormal- 
ization" prescription of Goulianos |2^ , which arises due to a breakdown of the triple- Regge 



theory for soft diffractive excitation. Also, nonfactorization is found in pre-QCD analysis 
of diffractive processes |2^ . An empirical analysis by Alvero, Collins and Whitmore pOf of 



the preliminary CDF double diffractive dijet data [§n|,0 indicates that hard factorization 
is violated in this process. In fact, their analysis suggests for parton distribution functions 
that are most consistent over all diffractive processes, the experimental DPE dijet cross 
section is much less than expected by factorization. On the other hand, the nonfactorizing 
CFS mechanism should enhance the cross section. Nevertheless, the analysis in ^^ does 



not rule-out experimental realization of the CFS mechanism, since general understanding 
from Regge models suggests that there is a large source of suppression which will emerge 
from effects generically termed absorptive corrections. These effect are due to exchanges 
of pomerons and gluons between particles in the basic model that possess very different 
rapidities, thus in particular between the two incoming hadrons. As such, these effects also 
are nonfactorizing. Actual computation of absorptive corrections is nontrivial since they 
are nonperturbative. Some work has been done to estimate their effects [ P^p3[ . A general 
conclusion of these works is that absorptive correction effects are independent of the hard 
kinematics and weakly s-dependent. As such, these effects should be very easy to distinguish 
from the N(L)DPE process. Also, these effects only should shift, in particular decrease, the 
values of the cross sections from those computed in our basic models and the effect should 
be the same for either the F(IS)DPE or N(L)DPE processes. In this paper, we are interested 
in examining qualitative differences between the F(IS)DPE and N(L)DPE processes, which 
are minimally model dependent. For this we will examine the Et and mean rapidity (?/+) 
spectra for both processes and for the cuts of CDF, D0 and representative cuts for LHC. 
We also will present total cross sections for all the models and all the experimental cuts. 
Thus the interested reader can test any suppression factor from any absorptive correction 
model that they wish. 

The reason that we do not give a demonstrative example of the overall absorptive cor- 
rection suppression factor is that, as will be seen in the sequal, for all the models that are 
examined, we find disagreement in the ratios of the cross sections from those found in the 
available experimental data. This discrepancy minimally is of order ~ 5. Present under- 
standing about absorptive corrections can not explain this discrepancy, since their effect 
only is to shift the cross sections by the same overall correction factor which drops out in 
the ratios. This discrepancy may reflect upon a limitation of our partonic level calculations 
or other controllable theoretical sources, or it may be that since the experimental data is 
still preliminary, it may yet be modifled. We will not attempt to formulate any theoretical 
explanations for this discrepancy found in this paper. Our modest goal is to examine the 
predictions of the basic models, which up to now still have not adequately been done for 



^Hereafter our usage of factorization without further specification always means in the context of 
hard factorization. 



these processes. Of special interest is to identify features that are minimaUy model depen- 
dent. Furthermore, in light of the breakdown of hard factorization suggested in ^|, it is 
important to know whether any features of the basic models are seem in the data. 

The paper is organized as follows. Sect. |I| reviews the kinematics of DPE dijet pro- 
duction and then models are presented for the nonfactorized and factorized processes. The 
inclusive dijet cross section also is reviewed in Sect. || and will be computed in later sections 
for comparison purposes. In Sect. |ITT| the CDF and D0 cuts for DPE dijet production are 
reviewed and representatives cuts for LHC are presented. In Sect. |^ results of our calcula- 
tions are presented for DPE and inclusive dijet Et and mean rapidity spectra and total cross 
sections. Subsect. |1V A| gives a general presentation of the results, Subsect. |IV B| examines 
the results in greater detail and Subsect. [IV Q compares our results with the preliminary 
data of CDF and D0 . Finally Sect. |V] presents our conclusions. Also in the last part of Sect. 
0, we discuss limitations of our models and compare with related models. Subsects. [IV B 



and [IV C present our results with considerable detail. For readers not wishing this much 
detail, the first part of Sect. concisely summarizes the basic results before proceeding to 
give our conclusions about them. 



II. MODELS 



This section reviews the kinematics of DPE dijet production and the formulas for the 



F(IS)DPE and N(L)DPE models, based on the presentation in |22|. 



A. Kinematics 



The DPE dijet process examined in this paper is shown in Fig. |l], 

p + p~^p' + p' + 2 jets. 



(1) 



The proton and antiproton collide at high center of mass energy s = {Pp + PpY —>■ oo, lose 
tiny fractions xjp/p and xjp/p of their respective longitudinal momenta, and acquire transverse 
momenta Qp and Qp. (This defines a diffractive regime, and in Regge theory would lead 
to an expectation of the dominance of double pomeron exchange (DPE)). Using light-cone 
coordinates (+, — ; _L), the components of momenta of the hadrons in Fig. |1] are 
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Here we use bold-face type to indicate two-dimensional transverse momentum. 



(2) 



The jets carry large momenta of magnitude Et in the plane perpendicular to the collision 
axis with azimuthal angle 0. (This defines a hard-scattering regime.) The small transfer 
of longitudinal momentum to the hard process implies large rapidity gaps between the jets 
and the two outgoing hadrons. The momentum delivered by the two incoming partons to 
the hard collision that creates the jets is some portion Xp^Xp of the longitudinal momentum 
fractions xjp/p, xjp/p respectively, < Xp < xjp/p, < Xp < xjp/p. Thus for the jets, ignoring 
terms of relative order -C 1, the components of their momenta are 



pi = iaxp J -,bxp J -; Et cos (f), Et sini 

P2= (bxp J -,axp J -^-Et cos (f),-ET sin (f) j , (3) 



where it is convenient to define 



l + Vl^^ 
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with 
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b=l-a, (4) 



K = . (5) 

JU-rjJU'pij 

For later use, define the ratios 

^ ^ ^^, ^ ^ ^^. (6) 

Xjpjp Xjpip 

It is conventional to describe the jet kinematics through the transverse momentum Et 
in Eq. (H) and the rapidity variables 

1 T) 1 T)~ 

which sometimes are expressed as y+ = (yi + y2)/2 and y^ = yi — y2. In terms of the jet 
rapidity variables and Et, we have 

xp = ^{e-y^+e-y^), (8) 

and 

f^ = 97 TT- (9) 

cosh2(y_/2) ^ ^ 



B. Factorized (Ingleman-Schlein) DPE - F(IS)DPE 

To obtain the expression for the Factorized (Ingelman-Schlein) DPE (F(IS)DPE) dijet 
differential cross section, first recall the inclusive dijet differential cross section (Fig. |]) 

7 dijet 

(J fT . 77" 

^ ^'"^^'^2 = 7 II[/i/p(%)/i/p(a^p)^ii(^> i^ ^) + fjlpMfi/p{xp)aij{s, i, m)]/(1 + 6ij), (10) 



^2 ., .. 



where s is the CM energy between the two incoming hadrons (here protons), fi/hi^h)-, is the 
inclusive parton distribution functions for parton species i in hadron h, and Xp, Xp are the 
parton momentum fractions relative to proton and antiproton respectively, carried by the 
two partons going into the hard interaction, dij is the parton 2 to 2 cross section for parton 



species i and j with explicit expressions given in ||34|. 

In the IS picture, they regard the pomeron as a hadronic particle. The pomeron is 
hypothesized to be created from the incoming proton and carries some momentum fraction 
Xip/hi h = p,p, of that proton's longitudinal momentum. In DPE hard expressions, one 
simply thinks of the collision of two pomerons in the same way as any two incoming hadronic 
particles. As such, the inclusive dijet expression, Eq. (|TUp above, applies to this case with 
two modifications. First s now must be replaced by the appropriate CM energy for the two 
pomerons, which is precisely xjp/pXjp/pS, where s here is the CM energy between the two 
incoming protons. Second a "pomeron flux factor" must be introduced, that expresses the 
probability to find a pomeron inside the proton. 

With these considerations in mind, the expression to F(IS)DPE dijet differential cross 
section is (Fig. H) 

I dijet 

— dxjp/pdxjp/pfjp/p{xjp/pjfjp/p{xjp/p 



dyidy2dE^ J "" ^/p—^/p^^/py ^/p^^^/pv '^"'' xjp/pXjp/pS 
^[fiM(^p)fj/Al3p)^ij{s, t, u) 

+ fJ/lp{Pp)f^MPp)^^A^,i,^)]/i^ + '^.j)- (11) 

In this expression Xp, Xp again are the momentum fractions of the incoming partons relative 
to the respective protons and (3h are the parton momentum fraction with respect to the 
pomeron, as defined in Eq. (||). fi/ip{Ph) now is the pomeron parton distribution function. 
dij is the parton 2 to 2 cross section, which is the same as in the above inclusive case Eq. 
(|10|). Finally, fip/h{xip/h) is the "pomeron flux factor" . In our work, we will use the pomeron 
flux factor of Donnachie and LandshoflQ |P5|| , 



^There is another commonly used pomeron flux factor which is of Ingelman and Schlein [^. 
This differs from the DL flux factor primarily in its normalization. However a change in the 
normalization factor completely is compensated for by changing the parton densities by an inverse 



factor. Thus the parton densities are obtained, for example in |36|, for a set of data without any 
a priori expectations as to their normahzation. 



f^u^lp) 



f^U^ip) 






Am; - 2.8t 
Arrip — t 



l-t/0.7^ 



X, 



l~2a{t) 
IP ■ 



(12) 



where nip ~ 0.938GeV is the proton mass, Po ~ l.SGeV is the pomeron-quark couphng 
and a{t) = ajp + 0.25t is the pomeron trajectory, ap is known as the pomeron intercept 



which for the soft Pomeron is ajp ~ 1.08 |37|. The pomeron parton densities used here 
are those of Alvero, Colhns, Terron and Whitmore (hereafter ACTW) ||3^. Their fits were 
to diffractive deep inelastic and diffractive photoproduction of jets, in which ajp was a free 
parameter that was fit to data and found to be ajp ~ 1.14. 

The ACTW fits are to five models, which covers a very general set of possibilities. 
Retaining their notation, the models will be denoted as ACTW A, B, C, D, and SG. The 
precise description of these models can be found in Sect. IID. of their paper. In brief, 
the models A-D use conventional shapes for the initial distributions. Model A represents a 
conventional hard quark parametrization, B has in addition to A an initial gluon distribution, 
C has in addition to A a soft quark distribution, and D has both additions to A. The final 
model, SG, has a gluon distribution that is peaked near (3 = 1. This form was motivated 
by the fit obtained by the HI collaboration. In [^, they refer to it also as the "super-hard 
gluon" . 



C. Nonfactorized (Lossless) DPE - N(L)DPE 

Our expression for the N(L)DPE dijet cross section is based on the toy quantum field 
theory model in ||22| which in effect is the model of Low-Nussinov-Gunion-Soper |lT7| , p!8 . 
The N(L)DPE dijet cross section expression obtained here extends from [^ to account for 
the one-loop Sudakov suppression factor. We presented preliminary results with Sudakov 



suppression in [^ 
and KhozeQ g 



Our treatment of Sudakov suppression is the same as by Martin, Ryskin, 
In fact, at one-loop order the nonabelian expression required here is the 
same as the abelian expression of Sudakov [^] which in the context of hard scattering was 
obtained earlier by Collins [^. The only difference is, the abelian expression must be 
multiplied by an overall group theory factor to account for the additional color degrees in 
the nonabelian case. 

For the N(L)DPE model in Fig. ^d, xip/p and xp/p again are the longitudinal momentum 
fractions lost by proton and antiproton respectively. In difference to the F(IS)DPE case, the 
momentum fractions for the incoming partons to the hard process are equal to those lost 
by the protons, Xp = xjp/p and Xp = xjp/p or equivalently Pp = (3p = 1. Qualitatively this 



^In 1^^ two types of double diffractive dijet expressions are given, which they call exclusive 
and inclusive. Both these expressions are nonfactorizing processes of the CFS type, with the 
exclusive case the same as our N(L)DPE model. Their inclusive case implements the same two 
gluon nonfactorizing mechanism as their exclusive case. The difference is, for the inclusive case the 
incoming protons can diffract to any final state, provided there are rapidity gaps between these 
final states and the hard process. This process is not relevant to this paper. 



means all the inomentuin lost by the diffractive protons is transferred into the hard process. 
This kinematics is similar to the superhard component reported for the case of single-sided 
diffractive dijet production by the UA8 pO|] . 

Our expression for the N(L)DPE dijet differential cross section is 



7 dijet 
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where 
M = 






gp(k, -Qp) gp{k, Qp) e,(k-Q,) e,(-k-Qp) A{i,j; ^Fsik', E^). (14) 



Here, the "polarization" vectors are defined as 



e.(k) 
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(15) 



cjh^k-, Q) are hadronic form factors with the explicit expressions of our model in Eqs. (10) 
- (12) of [^, and A{i,j;f) is the hard amplitude. Two hard subprocesses are possible, 
fi'i'5'2' -^ 9192 and 9i'92' -^ <li<i2- The calculations in [pT] , p2[ showed that the latter process 
gives zero contribution to the N(L)DPE dijet cross section when the final state transverse 
momentum of the two diffractive hadrons is zero. This should suppress quark jet production 
relative to gluon jet production. In |2^ this expectation was explicitly confirmed. Thus 
only the 91192' —>■ 9i92 hard process is relevant for N(L)DPE production.. The explicit 
expressions for A{i,j;f) for this process are given in Appendix A of p2[. The Sudakov 
suppression factor is Fsiii^yE^), which at one-loop order is [P],P^] 



Fs{k\E^) =exp -^(k- 



where 



^(k^E|) 



3«5CE|)j^2 

in 



2 7T.2 

Ulrp 



Erp 



4(k2 + /i2 



(16) 



(17) 



and /i (~ IGeV) is a low energy cutoff scale. 

The amplitude M. in Eq. (|14D is fixed up to an overall normalization which implicitly 
is specified through the hadronic form factors gn- Based on the same quantum field theory 
model, an expression for the hadron-hadron elastic scattering amplitude can be determined 
and that expression involves the same hadronic form factors cjh- In fact this expression for 



the elastic scattering amplitude essentially is the one of Low-Nussinov |T^ and Gunion-Soper 
||18| . Thus the free parameter in our N(L)DPE model that fixes the overall normalization 
of M. is chosen to yield the experimental value of the pp forward elastic cross section. The 
details of this procedure are given in 



III. EXPERIMENTAL CUTS 

This section reviews details about the CDF and D0 experiments that are relevant to 
the DPE dijet process. In particular, we state the cuts we will use to represent the CDF 
and D0 DPE dijet experiments^. Also, cuts are given that are representative cases for DPE 
dijets at LHC. 



CDF has presented results on double diffractive dijet production pO Jpll , ^^ at y^ = 
ISOOGeV with transverse jet energies Et > 7GeV. The experiment has one Roman pot 
on the — rapidity side, which detects the diffractive hadron, here p, going in this direction 
with the cuts 0.04 < xp/p < 0.095. On the + rapidity side, there is no Roman pot, only 
a rapidity gap requirement. Thus, in principle, there is no specified cuts on the outgoing 
diffractive hadron, here p, that goes in the + rapidity side. However, based on the rapidity 
gap length on this side, they obtain the estimate 0.015 < xjp/p < 0.035. The experiment 
places no explicit cuts on the dijet rapidity region. In our calculations, as our cuts, we will 
use the entire central detector region —4.2 < yi,y2 < 2.0. 

The CDF double diffractive dijet experiment has three shortcomings in its interpretation 
as the DPE dijet process. First, the lack of a Roman pot on the + rapidity side to detect the 
proton is a primary source of ambiguity in differentiating rapidity gap events that involve 
diffractive excitation of the proton versus pomeron exchange. Second, a heuristic guide for 
pomeron exchange is that in the diffractive event xjp/p^p < 0.05. Above this limit other 
Regge exchanges may be important or the interpretation as a diffractive process altogether 
may be questioned. Based on this guide, the CDF Roman pot cuts on the diffractive 
antiproton, 0.04 < xjp/p < 0.095, exceed the optimal region for interpretation as pomeron 
exchange. Third, CDF does not implement on-line jet triggering. Instead, they collect a 
sample of events that have a p track in the Roman pot and a rapidity gap in the +rapidity, 
p, side. From these events, they separate those cases in which there are two jets in the 
central region with Ej^ > 7GeV and discard the remaining data. Since generally jets are 
difficult objects to create and the central region typically is soft, most of the collected 
data is discarded. Furthermore, Ej- > 7GeV is a relatively low transverse jet momentum 
requirement for Tevatron jets and at this level jets are fairly cloudy objects that may be 
difficult to reconstruct and measure precisely. 

D0 has been examining double diffractive dijet production at two center-of-mass en- 



ergies ^/s = 630GeV and ISOOGeV [p5l -|48[|. Hereafter these two cases will be referred to 
respectively as D063O and D018OO. At present, the D0 experiment has no Roman pot^, 
so that DPE dijet production operationally is defined as two hard jets in the final state 
which are separated from both sides of the beamline by large rapidity gaps. The dijet cuts 
are Et > 12GeV for D063O and Et > 15GeV for D018OO with \yi\, jyal < 1-0 for both 



^Some years back the UAl also had reported on jet events with double rapidity gaps in pp-collisions 



at \/s = 630GeV |42|. However their reported results are insufficient to include in our analysis. 



^In Run II, which is expected to start in late 2000, D0 plans to have Roman pots on both sides 
of the beamline to detect both the proton and antiproton 143] 
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center-of-mass energies. 

The D0 approach has both an advantage and a disadvantage to the CDF approach. 
The advantage is D0 implements on-hne jet triggering. Thus, they are able to make a 
more efficient usage of their collected data sample. In addition 12 and 15 GeV jets are 
much better defined for identification by cone algorithms. The disadvantage of the D0 
approach is without Roman pots the same ambiguities experienced by CDF arise here and 
are magnified. In particular uncertainty remains about what portion of the double gap 
events are double diffractive, single diffractive/single pomeron exchange or double pomeron 
exchange. Moreover the momentum fractions lost by the proton and antiproton are not 
directly measured, so that in principle the experiment places no restriction on them. An 
upper bound on the momentum fraction lost by the p and p can be estimated by examining 
the maximum energy deposition in the hard event. From this, it can be inferred EH] that 



Xip/p,xip/p ~ 0.05. This bound is consistent with the heuristic notion that a pomeron can 
carry no more than ~ 5% of the diffractive hadron's momentum. For the D0 cuts used 
in our calculations, we will use the upper bound xjp/p,xjp/p < x^^^ = 0.05. The jet Et 
requirements imply a minimum energy must be deposited in the hard interaction region, 
which places a lower bound bound xp/p, xp/p > x^™ given by 



™min™max„ 



Er = f J ■ (18) 

For convenience, to accommodate both center-of-mass energy cases, we set xfi" = 0.001, 
since this limit is lower than the values given by Eq. (|TB|). 

For DPE dijets at LHC with ^/s = 14, 000 GeV, we use the following cuts. The transverse 
jet energy will be Et > lOGeV with the rapidity region —1.0 < yi,y2 < 1-0. These cuts 
represent standard expectations for jets in the hard DPE process. For the momentum 
fractions lost by the proton and antiproton, four cases will be considered, LHC-1: 0.002 < 
xiP/p,p < 0.03, LHC-1': 0.00006 < x^/p^p < 0.03, LHC-2: 0.002 < xip/p^p < 0.01, and LHC-2': 
0.00006 < xjp/p^p < 0.01. These cuts are estimates of where diffraction should be important. 
The upper bounds are slightly more conservative than the heuristic limit of 0.05. For the 
LHC-2 and LHC-2' cuts, where the upper bounds are 0.01, diffraction clearly should be 
dominant as supported by the ZEUS [0] and HI |]51| experiments. The lower bounds on 



xjp/p^p for LHC-1 and LHC-2 again are suggested by the Zeus and HI data. The lower bound 
on xjp/p^p for the LHC-1' and LHC-2' cases is based on the minimal energy condition for the 
hard interaction Eq. ([T^). Here, the lower limit of xjp/p^p > 0.00006 accommodates both 
cases, LHC-1' and LHC-2', since this bound is below x]p™ given by Eq. ([18|) . 

A. Summary of the Experimental Cuts 

For convenience, below the cuts we use to represent the various experiments are summa- 
rized. 

CDF: v^ = ISOOGeV, Et > 7GeV, -4.2 < yi,y2 < 2.0, 0.015 < x^/p < 0.035 (+rapidity 
side), 0.04 < xip/p < 0.095 (— rapidity side). 
D018OO: Vi = ISOOGeV, Et > 15GeV, -1.0 < yi,y2 < 1.0, 0.001 < xp/p^xp/p < 0.05 
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D063O: v^ = 630GeV, Et > 12GeV, -1.0 < yi,y2 < 1.0, 0.001 < xip/p,xip/p < 0.05 
LHC-1: ^ = 14,000GeV, Et > lOGeV, -1.0 < yi,y2 < 1.0, 0.002 < xip/p,xip/p < 0.03 
LHC-1': v^ = 14, OOOGeV, Et > lOGeV, -1.0 < yi, 2/2 < 1-0, 0.00006 < xip/p, xp/p < 0.03 
LHC-2: v^ = 14, OOOGeV, Et > lOGeV, -1.0 < yi,y2 < 1.0, 0.002 < xip/p,xip/p < 0.01 
LHC-2': v^ = 14, OOOGeV, Et > lOGeV, -1.0 < yi, 2/2 < 1-0, 0.00006 < x^/p, xp/p < 0.01 

IV. RESULTS 

In this section, the results of our calculations are presented. Then, various cross-checks 
and features of the results are discussed. Finally a qualitative comparison is made of our 
results with the preliminary results from the CDF and D0 experiments. 

A. Presentation 

Calculations have been performed of the ?/+-spectra, Et spectra, and total cross sections 
for the N(L)DPE and F(IS)DPE dijet processes with the CDF, D0 and LHC cuts that 
were discussed in Sect. |T|. For comparison, the standard inclusive dijet process, Eq. 
(p!0|), also has been computed with cuts comparable to the corresponding CDF, D0 and 
LHC DPE cuts. In particular, the corresponding inclusive dijet cuts we use are for CDF: 
v^ = ISOOGeV, Et > 7GeV, -4.2 < yi,y2 < 2.0, D018OO: y^ = ISOOGeV, Et > 15GeV, 
-1.0 < yi,y2 < 1.0, D063O: ^ = 630GeV, Et > 12GeV, -1.0 < yi,y2 < 1.0, and LHC: 
^/s = 14, OOOGeV, Et > lOGeV, -1.0 < |/i, 2/2 < 1-0. 

Our results for the y+ and Et spectra are presented in Figs. ^ - || for respectively CDF, 
D018OO, D063O and LHC. In all four of these figures, the (a) graphs contains the ?/+-spectra 
for the N(L)DPE and F(IS)DPE models, the (b) graphs contains the ratio of ?/+ spectra 
between the DPE and inclusive processes, {dacpE/ dy+) / {do'inci./ dy+) , the (c) graphs contain 
the i?T-spectra for the N(L)DPE and F(IS)DPE models and the (d) graphs contain the ratio 
of the Et spectra between the DPE and inclusive processes, {daopE / dET) / {damcij dET) ■ 
In the CDF and D0 figures, the solid curves represent the F(IS)DPE ACTW A-SG models 
and the dashed curves represent the N(L)DPE model with Sudakov suppression factor Eq. 
(p!6|) none (i.e. Fs = 1) and fx^ = 1.0, 0.3 GeV^. The D0 figures also have dashed-dotted 
and dotted curves, which represent for the F(IS)DPE and N(L)DPE processes respectively 
some modified cuts. The specifics of these curves will be explained at the appropriate time 
in the discussion that follows. For the LHC cuts in Fig. (g), the F(IS)DPE ACTW D model 
is represented for LHC-1, 2 by the solid curves and LHC-1', 2' by the dashed-dotted curves 
and the N(L)DPE Sudakov suppressed /i^ = 0.3 model is represented for LHC-1, 2 by the 
dashed curves and LHC-1', 2' by the dotted curves. 

The total cross sections for all the DPE dijet cases are in Table 1. For the corresponding 
inclusive dijet cases, the total cross sections are 

a^^f = 1.9 mb, 
^001800 = 0.023 mb. 



12 



al;^p = 0.71 mb. 

(19) 

For the inclusive process, the CTEQ5 parton distribution functions were used [52|. For 



the F(IS)DPE case, we use the best fit value of the pomeron intercept found by ACTW [^ 
ajp = 1.144. Also, the ACTW pomeron parton distributions were dependent on ajp as a 
consequence of their fitting procedure, and we have used the ones at ajp = 1.144. Both types 
of distribution functions are evolved with three fiavors of quarks, and in all calculations, we 
set Aqcd = 0.271GeV. 

Further details about the ACTW pomeron parton distribution functions can be found in 
their paper |3^. However some relevant facts about the ACTW results are reviewed here in 



order to put our calculations in perspective with their results. Amongst the five models of 
pomeron parton distribution functions considered by ACTW, they found that the ones with 
high gluon content gave the best fit to the diffractive DIS and diffractive photoproduction 
data. In their notation the high gluon models are B,D and SG, with D giving the best 
fit, whereas the low gluon or quark dominated models are A and C. As a cross check, the 



predictions of their fitted models were examined in ||29| for charm production in e-p collisions. 
The cross sections for the high gluon models were within an order of magnitude of both the 
ZEUS and HI data, with model D again doing the best, whereas the cross sections predicted 
from the low gluon models, A and C, were two orders of magnitude below the experimental 
data. Thus, hard factorization for diffractive lepton-hadron scattering appears to be well 
supported by the ACTW analysis. However, in confronting their fitted models to diffractive 
hadron-hadron scattering, a pronounced inconsistency with data occurs. For DPE dijet 



production, the analysis in [^ found that the high gluon model cross sections were 20-300 



times larger than the CDF data, with model D having the largest discrepancy, whereas the 

low gluon models actually agreed within a factor two of data. Other tests made by ACTW in 

3^ also revealed similar inconsistencies for hard factorization in diffractive hadron-hadron 



processes. 

At the moment there is no explanation for this breakdown in hard factorization, and 
before any insight may be gained, it appears the situation still is in search of more tests 
of the data. This paper provides several comparison tests between theory and experiment 
with the primary aim to discriminate between the two hard DPE mechanisms. However, en 
route, these tests also supplement the ACTW hard factorization analysis and may provide 
additional insight into the problems uncovered in their work. In particular, the dijet dis- 
tributions calculated in this paper provide more detailed predictions from the basic models 
than just total cross sections, with which to confront data. Furthermore, total cross sections 
are calculated for several experimental cuts, so that ratios amongst them can be tested to 
data. These ratios can help test the validity of the absorptive correction models, which, 
as discussed in the Introduction, in generic Regge physics inspired models are believed to 
be fairly independent of the hard kinematics. Thus, for example, for the CDF and D0- 
1800 cases, the absorptive correction effects from these models would give the same overall 
correction factor and so should cancel out in the ratio between the two experimental cross 
sections. 
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B. Discussion and Cross-Checks 

This subsection highlights some interesting features in the results and explains their 
underlying origins. The F(IS)DPE and inclusive cross section formulas for dijet production 
are well known in the literature. We simply will quote where necessary properties about the 
various quantities that enter in these expressions such as the parton distribution functions, 
pomeron flux-factor and hard matrix element. On the other hand, the N(L)DPE dijet cross 



section formula is less familiar. In [22] it was noted that for forward scattering of both 



hadrons and with no Sudakov suppression factor, Fs = I, the square of the amplitude Eq. 
(p!^) becomes 

|AT(0, 0)p = 647r (^] S,,6kiH,j,i, (20) 

where Hnjj is the square of the hard parton amplitude. Its exact expression is given in 



H, which evaluates to be Eij Hujj = 18[4:7ra{E^)]'^{N^ - 1) cosh^(?/_/2). Using these 



expressions, Eq. (|T3|) crudely can be approximated by the following expression which can be 
evaluated upon inspection. 



NimPE^iaa^ f rf^Q,rf^Q,- ^_^^'^":t ... 1^,(0,^1)^, (21) 



where the Sudakov suppression factor is approximated at the k = point. In this expres- 
sion, the integral of the two diffractive protons' outgoing transverse momentum phase space 
(Qp, Qp) can be approximated as 



r 



^2n ^2p)_ ^ 2|y.|2 /99\ 

'approximate ^ ^£^pUj ^^p '"^^ '' | '^ | yyiax ' \ / 

1— reeion 



where \t\max is a fixed parameter that represents the characteristic transverse momentum cut- 
off for the diffractive protons. With these approximations, and {da/dt)ei = 201mb/GeV^, 
which is obtained from the optical theorem from the total cross section afot = 62mb [|37| , 
Eq. (pi]) becomes 

7 dijet 2 / 7-i2 \ 

g|£g.(45.8.VGeV=)lC.^. (23) 

To determine |t|max, the above expression can be compared to the exact numerical expression, 
Eq. ([13D , at one point. From this, we will set \t\max ~ 0.26GeV^. 

The three subsections to follow examine the CDF, D0 and LHC cases in turn, with 
cross checks and explanations offered for the various features of the results in Figs. |^ - |^. 
One immediate cross check of our results is the magnitudes of the cross sections. For the 
F(IS)DPE case, we have verified that our results agree with [^. For the N(L)DPE caseQ, 



^The exclusive double diffractive model in [39| has a more detailed description of the two-gluon 
pomeron process compared to our model. For this reason, direct comparison of total cross sections 
is not possible between out model and theirs. In the Conclusion, we will discuss further the model 
inH 
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it will be seen below that the exact numerical results are consistent with the approximate 
expression Eq. (p3|). 



1. CDF 

The CDF results in Fig. ^ have the following noteworthy features. From Figs. |^a and 
^, the |/+-spectra for the N(L)DPE process (dashed curves) are localized to the region 
— 1 ~ I/+ ~ 0. The |/_|_-spectra is much broader for the F(IS)DPE (solid curves) versus 
N(L)DPE case, and from Fig. ^d both are less broad than the inclusive |/+-spectra. This 
difference in the broadness of the ?/_|_-spectra between the F(IS)DPE and N(L)DPE processes 
with CDF cuts is one of the most pronounced signatures found in this study that could help 
to differentiate the two processes. As will be seen below, this difference reflects upon intrinsic 
kinematic differences between the two processes, and thus is a reasonably model independent 
feature. Turning to the -E-r-spectra, from Fig. ^ the N(L)DPE process falls much slower 
than the F(IS)DPE process. In fact from Fig. ^d, the N(L)DPE process is seen to be 
almost fiat for Et < 45GeV for the two cases with Sudakov suppression and slightly rising 
for the case with no Sudakov suppression. Then for Et > 45GeV, all three cases rapidly 
fall to zero. For reference, an exactly fiat spectra in Fig. ^ would imply it has the same 
shape as the inclusive i^^-spectra. Thus the two Sudakov suppressed N(L)DPE E'-p-spectra 
have approximately the same shape as the inclusive E'-p-spectra. On the other hand, for the 
F(IS)DPE spectra, all five cases fall much more rapidly in Fig. |^d relative to the inclusive 
spectra, with the SG case falling the least rapidly. 

The behavior of the CDF N(L)DPE spectra can be understood from the approximate 
Eq. ( |25| ) and by examining the dijet rapidity phase space. Recall for N(L)DPE processes, 
the parton momentum fractions Xh that enter the hard interaction equal the corresponding 
pomeron momentum fractions Xjp/h, Ph = Xh/xip/h = 1- As such from Eq. (|^), the cuts 
on xip/h imply direct restrictions on the jet rapidities yi, ^2- One finds upon inspection of 
Eq. (H) and the explicit CDF rapidity cuts from Sect. |T| that at Et = 7GeV, dijets only 
appear in the rapidity ranges 1.4 ~ yi < 2, —3.3 ~ ^2 ~ —2.3 (and interchange yi < — > 1/2), 
which equivalently implies — 1.4 ~ y+ ~ —0.15. As Et increases, the kinematically allowed 
rapidity bands move inwards towards zero rapidity. Generally both bands also get narrower. 
However, since the rapidity band at the proton side (+rapidity) was prematurely cut-off at 
2 due to the explicit rapidity cuts, this band first broadens up to Et ~ 14GeV and then 
narrows thereafter for higher Et- As such at Et = 14(28) GeV dijets appear in the bands 
0.5(-0.3) ~ yi ~ 1.5(0.75), -2.6(-1.8) ~ ^2 ~ -1.5(-0.7) (and interchange yi < — > 2/2), 
which corresponds to — 1.1 ~ y_|_ ~ at both Et scales. These considerations suffice to 
explain the localized y+-spectra in Figs. §a and ^o for the N(L)DPE process. 

Applying these estimates to Eq. (p3D , the differential cross section daj\[(^QDPE/dET at, 
for example, Et = 7, 14, 28GeV gives the relative magnitudes 1, 0.2, 0.03, which are within 
a factor 2 of the exact numerical results in Fig. ^. This figure also indicates that the region 
above Et > 30GeV accounts for less than 0.5% of the total N(L)DPE cross section. A 
check of the dijet phase space indicates that above this Et, the accessible region is rapidly 
diminishing. In fact due to the kinematic constraints, the maximum energy that can be 
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deposited in the hard region for either the F(IS)DPE or N(L)DPE processes is Jx^%Xp%s 

which for the CDF cuts imphes the largest dijet Et is Et ~ 51GeV. For the N(L)DPE case, 
this cutoff is best seen in Fig. ^. 

The last point to address about the CDF N(L)DPE process is the magnitude of the 
total cross section. The exact numerical results are given in Table 1. Estimates based 
on Eq. (p3D, where the phase space integral and all other quantities are approximated at 
Et = -E™™ = 7GeV, agree up to a factor 2 with the results in Table 1, including the ratio 
amongst the three cases of Sudakov suppression, none, /i^ = IGeV^, and fi"^ = 0.3GeV^, of 
respectively 1, ~ 0.5, ~ 0.25. 

Turing to the F(IS)DPE process, the first point to be addressed is the steeper decrease 
of the Et- spectra in Figs. §c and ^i relative to both the N(L)DPE and inclusive E^- 
spectra. Two facts are useful for this analysis. First, as Et increases, in general, the 
average value of Xp and Xp increase, since more energy must be deposited into the hard 
region. Second, the pomeron parton distribution functions fi/ip{l3h) at small argument 
Ph = Xh/xp/h grow as fi/ip{(3 — > 0) ~ /5~" with ~ a ~ 1.5, and at large argument 
vanish as fi/ip{P -^ 1) ~ (1 — /?)'' with ~ 6 ~ 1.0. Thus at small Et, Xp and Xp, and 
so therefore also Xp/xjp/p and Xp/xjp/p, are closer to their kinematic lower bounds, which 
implies the parton distribution functions are at their largest. However as Et increases, it 
imphes Ph = Xh/xjp/h — ^^ 1 so that fi/ip{xh/xjp/h) — > 0. In contrast, within the same Et 
range, the behavior of the inclusive parton distribution functions is very different, primarily 
due to the difference in behavior of their arguments. The inclusive distribution is evaluated 
with respect to Xh not Xh/xjp/h- Since x™™ < Xh < xjp/h and xjp/h < 0.1, Xh within 
this range always is relatively small. Thus the inclusive parton distributions within the 
equivalent Et range have less variation and generally are large. This difference in behavior 
of the arguments for inclusive and F(IS)DPE parton densities explains the steeper decline 
of the latter's E'r-spectra relative to the former. 

There are two immediate checks that verify the above observations about the ii^r-spectra. 
First pomeron parton distribution functions that fall slower as /5 — ;> 1 should have flatter 
i?r-spectra in Fig. |^d and this is the case for the ACTW SG model. Second, if the upper 
limits on Xjp/pp are increased, then for fixed Xh, the ratio Xh/xjp/h is smaller. Thus, for the 
same jet kinematics, fi/ip{xh/x]p/h) should fall less rapidly, which in turn would flatten the 
ii^r-spectra in Fig. |^d. One can verify this effect for any general pomeron parton distribution 
function. 

To further quantify the above observations about the CDF case, we can ask how small 
do the the parton momentum fractions become. From Eq. (0), naturally the minimum 
value for one occurs, when the other is at its maximum. As a more realistic estimate, let us 
assume the "large" region for the parton momentum fractions is when Xh/xjp/h ~ 0.5, since 
above this point the pomeron parton distribution functions rapidly vanish. Thus the parton 
carrying the "large" momentum fraction will have Xh^^^ ^ ~ 0.5x^/^ . We substitute Xhi^^ ^ 
in the LHS of Eq. (pISl). We now ask above what Et will the other parton momentum fraction 
also be in the "large" region, under the assumption that when both parton's momentum 
fractions are "large", there is negligible contribution to the cross section. By this criteria, 
we find that both momentum fractions are in the 'large" region, which means Xp ~ 0.025 on 
the proton side and Xp ~ 0.03 on the antiproton side, once Et ~ 35GeV. 
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Although for the pomeron case, these momentum fractions are large because /3/i ~ 1, the 
situation is different for the inclusive case. The inclusive parton distribution functions are 
evaluated with respect to Xh, not Ph = Xh/xjp/h- As such the range for their arguments is 
Xp,P ~ 0.02 — 0.04 and within this range the parton distribution functions have very little 
variation. The numbers quoted in this example are very crude, but they illustrate the reason 
in Fig. ^ for the F(1S)DPE Et spectra's steeper decline relative to the inclusive case. 

The above discussion ignored entirely complications from the pomeron flux factor f^^. 
This is because its approximate behavior is f^^ ~ Xp "^ and in either of the two ranges 
0.015 < xp/p < 0.035 or 0.04 < xp/p < 0.095, its variation is relatively small, i.e. less than 
a factor 3. 

It is worth noting that the rise in the N(L)DPE Et spectra in Fig. |^d has similar 
explanation to the one given above for the differences in ^^^-spectra between the F(IS)DPE 
and inclusive cases. In particular, the inclusive parton distribution functions will fall a 
little as Ejn increases since the average values of Xh will increase. However the proton form 
factors in our N(L)DPE model are insensitive to this variation. This is one of the notable 
differences between the N(L)DPE model and both the inclusive and F(IS)DPE models, and 
it explains the relative rise in the former's Et spectra to the latter. Furthermore, the rise is 
less pronounced for the Sudakov suppressed N(L)DPE processes, since they provide greater 
suppression to the N(L)DPE differential cross section da/dET as Et rises. 

In the N(L)DPE model, this lack of x/j- dependence in the proton form factors is not a 
fundamental requirement for nonfactorization of the CFS type. This is a simplying limitation 



in this particular model. In ^9|] nonfactorizing models similar to our N(L)DPE model were 



treated, except with more detailed modeling of the Xh dependences. We will discuss the 



models in |^ later in the paper, but we will not explore such modification to our model in 
this paper. 

Next, we will understand the behavior of the j/+-spectra for the F(IS)DPE case in Fig. §a 
and ^. To simplify the problem, we assume the main features of the spectra are determined 
by the low Et dijets Et ~ _E™™ = 7GeV. We want to understand why the j/+-spectra is 
much broader for the F(IS)DPE versus N(L)DPE case and moreover why the former also 
is skewed towards the — j/+ side. For this, note that at fixed Xp, Xp, and Et, the largest y+ 
attainable by the dijets is when both are on the same side with equal rapidity, y_ = 0. In 
this case Eq. (||) becomes 

XP = ^exp{±y+). (24) 

By evaluating Xp^p within their allowed range, this expression gives the limits on y+. The 
allowed ranges for Xp^p are determined by the same criteria as before that both momentum 



fractions must be "small", Xh/xjp/h < 0.5. This condition implies the ranges Xp ~ 0.02, Xp ~ 



0.05, with the lower limits in both cases governed by the energetic condition Eq. ([T8|) . By 



these crude approximations, the y^ range for the spectra is —1.9 ~ ^+ ~ 1.0, which coincides 
reasonably well with the exact numerical results in Fig. ^. Furthermore, one finds as the 
boundaries of the j/+ range are approached, Xp and Xp are increasing for both the inclusive 
and F(IS)DPE cases. As such, the contributions from the respective parton distribution 
functions are decreasing. In particular, the distribution functions decrease slower for the 
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inclusive versus F(IS)DPE case, since the range of the argument in the former is much 
smaller ~ 0.1 versus the latter ~ 1. This part of the explanation is the same as our earlier 
discussion which compared the i?T-spectra for these two processes. The final outcome is the 
F(1S)DPE |/+-spectra for all the models in Fig. |^ are narrower than the inclusive one. 

The last point to note is that the total cross sections in Table I for the five models come 
in the ratio 1(A):100(B):1(C):400(D):20(SG). To obtain insight into these ratios, it is useful 
to decompose the cross section in terms of the parton initiated processes. For the B,D and 
SG models we find ~ 80% of the cross section comes from the pure gluon process dgg and the 
remaining fraction predominately from the dgq process. On the other hand, for the A and C 
models, the cross section decomposes as less than 5% from o"„„, ~ 40% from a„„ and ~ 60% 







from the pure quark initiated process (Tqq + dqq. Therefore the B,D, and SG models are more 
gluon controlled whereas the A and C models are more quark controlled. However, for none 
of the five models is it the case that one species of partons, quarks or gluons, dominates the 
cross section. 

To cross check these findings, we examine the pomeron parton distribution functions in 
the most probable [3 range, which we estimated above to be ~ 0.01 — 0.1. Two features 
about the parton distribution functions are evident. First, in this (3 range, the A and C or 
B and D distributions are the same magnitude, with the latter pair about a factor 10 larger 
than the former pair and the SG distribution is a factor 3-5 larger than the former. Second, 
for a given parton model, the ratio of the gluon to quark parton distribution function in this 
13 range, fg/ip{P)/ fq/ip{l3), is for the A and C models 3 — 20, the B and D models 20 — 50 
and the SG model 10 — 30. These properties are consistent with the general trends found 
amongst the five models for the total cross sections in Table 1. However since this turns 
out to be an intermediate regime between quarks and gluons, these simple indicators are 
insufficient to better quantify the results found from the exact numerical calculations. 

2. D0 

The D018OO results in Fig. ^have the following interesting features. Both the N(L)DPE 
(dashed curves) and F(IS)DPE (solid curves) |/+-spectra are localized to the region \y^\ < 1. 
In contrast to the CDF case, here the F(IS)DPE |/+-spectra is slightly narrower than the 
N(L)DPE |/+-spectra, which best is seen in Fig. |]b. For the ii^T-spectra, similar to the the 
CDF case the N(L)DPE process falls with increasing Et much slower than the F(1S)DPE 
process. In comparison to the inclusive process, for the N(L)DPE case both the ?/+-spectra 
in Fig. ^D and E^ spectra for Et ~ 40GeV in Fig. ^ are fiat, thus have the same shape 
as the corresponding inclusive spectra. On the other hand, for the F(IS)DPE case, the ?/+- 
spectra in Fig. |^ is much more localized than the inclusive spectra and the ii^7--spectra in 
Fig. pi falls much faster. 

The basic features of the N(L)DPE spectra can be understood, once again, through the 
approximate cross section formula Eq. (p3|) and by examining the available dijet rapidity 
region based on the explicit D0 rapidity cuts and those implied by the cuts on xjp/p^p. 
Carrying out this analysis, at the lowest E^ = 15GeV we find that the cuts on xjp/p^p 
place no additional restrictions, so that the available rapidity region is —1 < yi,y2 < 1- 
As Et increases, the first rapidity region that diminishes is for same-side jets and starting 



at the periphery. The reason is evident from Eq. (||). At fixed \yi\ and I1/2I, Xp and Xp 
differ more for same-side dijets than for opposite-side dijets. As such, the larger of the two 
parton momentum fractions will reach its upper limit at smaller yi and 1/2 for same-side 
versus opposite-side dijets. For example, the available rapidity region at Et = 20(30)GeV 
for opposite-side jets is unchanged < yi < 1 and — 1 < ^2 < (and interchange 
1/1 ^^ 1/2) whereas for same-side jets the rapidity regions are < yi,y2 < 0.8(0.4) and 
— 0.8(— 0.4) < yi,y2 < 0. By Et > 30GeV all regions of rapidity space diminish. For 
example at Et = 40GeV, the allowed rapidity space is < yi ~ 0.5, —0.5 ~ 1/2 < (and 
yi ^->- Z/2) and < 2/1,1/2 ~ 0.2 and —0.2 ~ 2/1,2/2 < 0. Finally for Et ~ 45GeV there is no 
allowed rapidity region. 

Integrating over the rapidity region in Eq. ( ^3]) with these estimates, we find consistency 
with the exact numerical results for da{ET)/dET in Figs, ^c and ^. Also the rapid cutoff in 
the i?T-spectra at Et ~ 45GeV, best seen in Fig. ^, is consistent with our crude estimates 
here that the available rapidity phase space vanishes at this point. 

For the N(L)DPE 2/+-spectra in Fig. |a, we again can apply the above results, except 
integrating Eq. (^) over Et and 2/_. The basic shape of the 2/+-spectra can be understood 
by assuming dominance of the low Et regime Et ~ 15GeV and examining the behavior 
of the 2/_ phase space as a function of 2/+. For the latter we find at 2/+ = the 2/- range 
is A2/_(2/+ = 0) = 4 and it vanishes linearly to zero as 2/+ — *■ 1- In fact relative to linear 
decrease as a function of 2/+, the 2/+-spectra in Figs, ^a and §3 is more enhanced near the 
middle 2/+ = relative to the periphery 2/+ = 1- This is an effect of higher i^^-dijets. Recall 
from above that as Et increases, the rapidity region first to diminish is for same-side dijets, 
or equivalently the large 2/+ region. For example at Et = 20GeV the region 1 < I2/+I < 0.8 
no longer has dijets. In general higher Et dijets enhance the region near 2/+ = relative to 

To estimate the N(L)DPE total cross from the approximate expression Eq. (p3D, the Et 
integral is performed with the rapidity phase space region evaluated at Et = 15GeV. The 
latter yields a rapidity area ~ 4, which implies for no Sudakov suppression the estimate 
a ~ 1.7/ib. Evaluating the Sudakov suppression factor at Et = 15GeV leads to suppression 
factors to the total cross section of ~ 0.67 and ~ 0.11 for /x^ = IGeV^ and /x^ = O.SGeV^ 
respectively. These estimates are within a factor 2 of the exact numerical results in Table 1. 

Turning to the F(IS)DPE process, the steeper ii^-r-spectra relative to both the N(L)DPE 
and inclusive processes arises for reasons similar to the CDF case discussed earlier. In par- 
ticular, the high Et region in general requires higher parton momentum fractions, which are 
suppressed by the pomeron parton distribution functions. For Et ~ 22GeV both incom- 
ing partons to the hard process carry "large" momentum fractions, Xh/xjp/h > 0.5. Thus 
above this Et, one should expect the Et spectra to diminish as evident in Figs. ||c and |d. 
For example relative to the inclusive i^T-spectra, the F(IS)DPE ii^^-spectra decreases for the 
A,B,C,D, SG models at Et = 22GeV by the additional factors 6, 9, 6, 9, 3 respectively and at 
Et = 30GeV by the additional factors 100, 250, 100, 250, 27 respectively. This faster decline 
for the F(IS)DPE Et spectra arises because the edge (i.e. /3 ~ 1) of the pomeron parton 
distribution functions is being reached. Consistently, the slowest decline of the ii^T-spectra 
amongst the five models is by the SG model, whose parton distribution function decreases 
the slowest at /3 — ;> 1. 
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For the ?/+-spectra in Figs, ^ja and ||b, their shape can be understood through rapid- 
ity phase space consideration and the behavior of the pomeron parton distribution func- 
tions. Generally larger \y_^\ requires larger Xp and Xp. Similar to the N(L)DPE case, at 
Et = 20GeV, for example, the xip/p^p and dijet rapidity cuts prohibit dijets for |y+| ~ 0.8. 
In addition to simple phase space restrictions, for the F(IS)DPE process larger parton mo- 
mentum fractions Xp and Xp thus larger \y+\ are further suppressed due to the pomeron 
parton distribution functions. This additional source of suppression at large \y^\ explains 
in Figs. ^ and ^d. why the F(IS)DPE |/+-spectra are a little narrower than the N(L)DPE 
ones. 

Turning to the magnitude of the cross sections, from Table 1 the relative sizes amongst 
the five parton distribution functions are about the same as in the CDF case discussed 
earlier, 1{A) : 50 (-B) : 1(C) : 200(D) : 50{SG). These ratios can be cross checked with the 
general behavior of the F(IS)DPE cross section formula, similar to our earlier treatment for 
the CDF case. Basically one finds that these ratios are consistent with the ratios of the 
pomeron parton distribution functions in the typical f3 regime (/? ~ 0.1). 

Finally, it is interesting to compare cross section magnitudes between D018OO and CDF. 
If the DPE dijet process is N(L)DPE dominated, the CDF and D018OO cases differ primarily 
by phase space area and the ^/E^ scaling factor. Despite the larger CDF central rapidity 
region, once their xp/p^p cuts are considered, the dijet phase space area between CDF and 
D018OO is approximately the same. As such the predominant difference between the CDF 
and D018OO cross sections is due to the l/E'^ factor, which if evaluated at their respective 
minima implies the N(L)DPE cross sections of D018OO should be a factor ~ 4 smaller than 
those of CDF. This is consistent with Table 1. 

On the other hand, for a F(IS)DPE dominated dijet process, in addition to the above 
two factors, an additional difference arises from the pomeron parton distribution functions. 
They make a significant difference due to their rapid growth at small-/?. Very roughly, since 
for CDF E"™™ is a factor 2 smaller and the average xjp is a factor 2 larger compared to 
D018OO, the typical j3 at which the pomeron parton distribution functions are evaluated in 
the CDF case should be a factor 4 smaller relative to D018OO. By knowing the typical (3 
range for the two cases, we can estimate the behavior of the parton distribution functions 
in that range. We expect that within the kinematically accessible range of /3, the typical 
range that dominates the cross section will be near its minimum limit since that is where 
the parton distribution functions will be the largest. The smallest (3 possible based on Eq. 
(0) is 0.015 and 0.11 for CDF and D018OO respectively. However at this limit, although 
one parton distribution function will be large, this effect is compensated by the other parton 
distribution function, which must be evaluated at /3 = 1 where it vanishes. Thus (3 regions 
away from this limit also will contribute significantly. As an estimate of an upper bound to 
the typical (3 range, we estimate at i?™™ for Xp = Xp and when xjp is at an average value, 
which we will take as ~ 0.05 and ~ 0.025 for CDF and D018OO respectively. Then we 
obtain Xh = 0.0077 and 0.017 so that (3 = 0.15 and 0.70 for CDF and D018OO respectively. 
So in summary the typical (3 ranges are 0.15 — 0.015 for CDF and 0.70 — 0.11 for D018OO. 
These estimates confirm that the typical f3 are a factor 4 — 7 smaller for the CDF case 
versus the D018OO case. Moreover the size of (3 for both cuts is O(10~^). In this (3 range, 
the growth of all five ACTW pomeron parton distribution functions is ~ l//5^'^~^'^. Thus 
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the typical size of a poineron parton distribution function in the CDF case will be a factor 
~ (4 — 7)i-2-i-3 i^igggj- than the D018OO case. Accounting for this factor for each of the two 
parton distribution functions and the factor 4 from Et scaling, we expect the cross section 
for any given pomeron parton distribution function model to be a factor ~ 10^ larger for 
CDF relative to D018OO. This essentially is what is found in Table 1. 

Note that the size of the typical /3 found above is interesting in its own right. It implies 
the range of f3 probed in both the CDF and D018OO cases is not tiny. Recall that very tiny 
(3 is the natural regime for gluon dominance. Thus, for the CDF and D018OO cuts, one 
should not assume that gluon dominance is necessary and we have shown earlier by explicit 
examples that there are parton distribution functions models (A and C in particular) in 
which that assumption is wrong. 

For D063O, the qualitative features of the spectra basically are the same as for D0- 
1800. In Fig. ^, the ?/+-spectra for the N(L)DPE process is a little broader compared to 
the F(IS)DPE process. For both the F(IS)DPE and N(L)DPE processes, the j/+-spectra is 
contained in a much smaller region, \y^\ ~ 0.3, compared to D018OO. However, similar to 
D018OO, the E'j'-spectra for the N(L)DPE case falls much slower than for the F(IS)DPE case. 
From Fig. |^, the relative decline of N(L)DPE E'r-spectra to the inclusive case is faster than 
in the D018OO case. More noticeably, in comparison to D018OO, the F(IS)DPE ii^^-spectra 
falls much faster relative to the corresponding inclusive -E^-spectra. In particular, the ratio 
of the F(IS)DPE to the inclusive ii^^-spectra for D063O falls by two orders of magnitude 
within an increase of Et by 2GeV whereas for the D018OO case the same decrease requires 
an increase of Et by 15GeV. 

To understand the features of the N(L)DPE case from Eq. (|23|) , first note that at 
Et = 12(15)GeV the allowed rapidity region for opposite-side dijets is < yi < 0.75(0.40), 
— 0.75(— 0.40) < ^2 < (and yi ^^ 7/2) and same-side dijets is < yi,y2 < 0.25(0.10), 
— 0.25(— 0.10) < yi,y2 < 0. By Et > 18GeV there is no accessible jet rapidity region. 
As such the explicit jet rapidity cuts of \yi\ < 1 are irrelevant since the upper bounds 
xip/p,p < 0.05 already prohibit sufficient energy deposition to produce dijets at the higher 
\yi\ region at even the lowest permissible Et- These crude estimates along with Eq. (p3D 
show consistency with the exact numerical results in Fig. |^ and Table 1. Finally since 
the jet rapidity region is rapidly shrinking even at the lowest Et, this explains the more 
rapid decline of of the E'r-spectra between the N(L)DPE and inclusive processes for this 
case compared to D018OO. This effect can be reversed by increasing the upper limit on 
xip/p^p- For example, the dotted curves in Figs. |^a-d present the various N(L)DPE spectra 
for Sudakov suppression with /i^ = 0.3 and when xjp/p,xjp/p < 0.1 

For the F(IS)DPE process, the explanation for the two spectra essentially is the same as 
in the D018OO case. The important quantitative difference is the typical parton momentum 
fractions are much bigger here than for D018OO. This explains the sharper fall of both the y+ 
and Et -spectra. For example, to produce 12GeV dijets with the minimal energy deposition, 
so at y_ = 0, requires for symmetric parton momentum fractions Xp = Xp = 0.038. Thus 

almost all the D063O events are in the "large" /? regime, Ph = Xh/xjp/h ~ 0.5. Recall in 
this region the pomeron parton distribution functions are rapidly diminishing. As such, 
the sharp decline in the i^T-spectra for D063O relative to the inclusive case primarily is 
because the periphery of the parton distribution functions at /3 ~ 1 are being probed. This 
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sharp decline can be reduced if the hmits on xjp/p and xp/p are increased. For example, the 
dashed-dotted curves in Figs. |^-d present the various F(IS)DPE spectra for the ACTW 
model D when Xjp/p, Xip/p < 0.1 

3. LHC 

The LHC results in Fig. |^ have the following interesting features. In Figs. |^ and 
Pb, the ?/+-spectra for the N(L)DPE LHC-1 and LHC-2 cases (dashed curves) are much 
narrower than for all the other cases. Relative to the inclusive |/+-spectra in Fig. Pd, all 
the F(IS)DPE cases and the N(L)DPE LHC-1' and LHC-2' cases (dotted curves) are flat, 
whereas the N(L)DPE LHC-1 and LHC-2 cases drop-off as |y+| increases with the latter 
falling fastest. For the Ex-spGciia. in Fig. ^jc, the most intersting feature is for the N(L)DPE 
LHC-1 and LHC-2 cases, for which the ii^^-spectra actually first rises with increasing Et 
until Et ^ 15GeV and thereafter falls. In contrast for N(L)DPE LHC-1' and LHC-2', the 
i?r-spectra are of a more standard behavior. For the F(IS)DPE i^^-spectra, the LHC-2 
(solid curves) and LHC-2' (dashed-dotted curves) cases fall much faster with increasing Et 
than the LHC-1 (solid curves) and LHC-1' (dashed-dotted curves) cases. In fact, in Fig. Ijd 
the F(IS)DPE LHC-1 and LHC-1' Et spectra are almost as fiat as the N(L)DPE LHC-1' 
and LHC-2' cases. 

For the N(L)DPE case, the primary difference between the primed and unprimed spectra 
arise due to the different lower bounds on xp/p^p of 0.00006 and 0.002 respectively. Due to the 
high a/s relative to the CDF and D0 cases studied earlier, much smaller parton momentum 
fractions are necessary to produce kinematically identical dijets. In fact the larger lower 
cut-off in xip/p^p for the LHC-1 and LHC-2 cases already is too large to produce substantial 
numbers of dijets in the range lOGeV < Et ~ 15GeV. For example, the only rapidity 
region where Et ~ lOGeV dijets can appear is for opposite-side jets with 0.6 < yi < 1, 
— 1.0 < 1/2 ~ —0.6 (and yi <-> 7/2). Once Et rises to 15GeV, the same-side jet region becomes 
accessible starting with the region yi ~ ^2 ~ and moving outward to higher rapidity with 
increasing Et- This behavior also explains the narrower ?/_|_-spectra for the N(L)DPE LHC-1 
and LHC-2 cases in Figs. ^ and ^. In the dominant Et range lOGeV < Et ~ 15GeV, 
dijets predominately emerge within ~ 0.2 ?/+-rapidity units about y+ = 0. The shoulder of 
the ?/+-spectra in Fig. ||a at t/+ ~ 0.4 corresponds to Et ~ 20GeV. 

In contrast, for the LHC-1' and LHC-2' cuts, the lower limit on xjp/p^p imposes no con- 
straints on jet rapidity. In this case, starting at Et = lOGeV, the complete rapidity region 
— 1 < 2/1, 2/2 < 1 is accessible. The explanation for the ?/+ and Et spectra in this case follows 
similar reasoning to the D018OO cases discussed earlier. 

To compare the magnitude of the N(L)DPE cross sections in Table 1 with Eq. (p3D, 
the dijet phase space must be estimated. For the N(L)DPE LHC-1' and LHC-2' cases at 
Et = Et™ = lOGeV the entire explicit rapidity region —1 < 2/1,2/2 < 1 is accessible with 
no additional constraints from the xjp cuts. This obtains the estimate with no Sudakov 
suppression of a ~ 4/ib. For the N(L)DPE LHC-1 and LHC-2 cases, because of the the 
large lower limit on xjp, at Et = lOGeV dijets only appear in a small region of opposite-side 
dijet rapidity space with Ay+Ay_ ^ 0.3. However at Et ~ 15GeV, the entire opposite-side 
jet rapidity region is accessible 0<2/i<— 1, — 1<2/2<0 (and 2/1 ^^ 2/2) but as yet only a 
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negligible region of same-side jet rapidity space. Thus the total rapidity area is Ay+ Ay_ ^ 2. 
Observe that the gain in phase space area is a factor 3 greater than the suppression factor of 
9/4 from the 1/ E"^ behavior of a. This is consistent with the rise in the i^y-spectra in Fig. 
||a. Applying the above estimates to Eq. (p3D implies for no Sudakov suppression a ~ 1/ib. 
The Sudakov suppression factor evaluated at Et = lOGeV implies the cross section in both 
cases should decrease by factors 0.4 for fj!^ = IGeV^ and 0.2 for yU^ = 0.3GeV^. These crude 
estimates are consistent with the exact numerical results in Table 1. 

For the ?/+ F(IS)DPE spectra in Fig. ^b, they are very similar to their inclusive coun- 
terparts. This contrasts the CDF and D0 cases in Figs. ^, ^b, and |^. The reason is the 
higher y/s, which for fixed Et and ?/+, requires smaller parton momentum fractions Xp^p. 
Assuming the shape of the |/+-spectra is dominated by the lowest Et ~ lOGeV region, for 
the full range of |y+| < 1, /5j = Xh/xjp/h typically never is "large" /?« < 0.5. Thus for the full 
I/+ range the pomeron parton distribution functions typically are not probed near /? — »• 1 
where they vanish. Instead, they typically are probed at intermediate jS regions, similar to 
the situation in the corresponding inclusive case. This is the basic reason for the similarity 
in the |/+-spectra between the F(IS)DPE and inclusive processes. 

For the F(IS)DPE Ey-spectra from Fig. |d the flatness of the LHC-1 and LHC-1' 
cases again arises from the small typical (3 values at which the pomeron parton distribution 
functions are evaluated. The sharper decrease of the F(IS)DPE LHC-2 and LHC-2' cases 
arises due to the lower upper bounds on xjp/p^p of 0.01, versus 0.03, for the LHC-2 and 
LHC-2' cases. The main consequence of these different upper bounds appears in the parton 
distribution functions and not from phase space. This is evident since only the latter effect is 
relevant for the N(L)DPE cases, and for them the E'r-spectra at large Et has no pronounced 
difference in all 4 cases. On the other hand, for the pomeron parton distribution functions, 
the two different upper bounds on Xjp/pp imply that for fixed Et, they are probed in the 
LHC-2 and LHC-2' cases at typical j3 values that roughly are 3 times larger relative to the 
LHC-1 and LHC-1' cases. This factor 3 difference has nonnegligible effect. Noting that the 
order of magnitude of the typical Xh is > O(10~^) or equivalently Ph > 0.06, with this lower 
bound increasing with Et, the factor 3 means the difference between (3 close to 0.1 where 
the parton distribution functions are sizable and /5 closer to 1 where they vanish rapidly. 

C. Comparison with Experiment 

Both CDF 0,§I],|3|,||,g and D0 [||-|§ have reported preliminary results on the 



double diffractive dijet process. In Sect. |T| limitations of both experiments were discussed 
which prohibit explicit measurement of the DPE dijet process in the present runs. Neverthe- 
less, one likely possibility considered by both experiments is that the majority of the double 
diffractive dijet events were DPE. Under this assumption, some of the reported features from 
these experiments will be interpreted below in terms of the models examined in this paper. 
D0 reports its i?T-spectra for the DPE dijet process to be similar to the inclusive Et- 



spectra with comparable cuts E5|-E3]. The earlier CDF preliminary reports |30J3l|H3H3 



also found this, although their most recent report |^ finds the ii^T-spectra falls much faster 
than for the comparable inclusive ii^r-spectra. Based on Figs. ^, ^J^d,^, the D0 and ear- 
lier CDF ii^T^-spectra lean towards that for the N(L)DPE process. Amongst the F(IS)DPE 
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processes, the E'r-spectra for the best fit ACTW model, D, as well as model B differ signif- 
icantly from their inclusive counterpart. However they could be consistent with the most 
recent CDF report [Q. The F(IS)DPE model with the most similar i^^^-spectra to the 



inclusive process, is ACTW SG. Recall in this model the gluon density peaks near P = 1. 

In the D0 case, as noted earlier, by increasing the upper bounds on X]p/p and xp/p, the Et 
spectra could be more flattened in Figs Ijd and 0d. In Figs. ^ and |^, the various spectra with 
upper bounds xjp/p^xjpjp < 0.1 are shown as the dotted curves for the N(L)DPE case with 
Sudakov suppression at /x^ = 0.3 and the dashed-dotted curves for the F(IS)DPE ACTW D 
model. Observe the dashed-dotted curves in Figs. ^ and [7|d for the ACTW D model are 
significantly flatter than any of the F(IS)DPE cases with upper bounds xjp/p^p < 0.05 (solid 
curves) (for comparison, the dotted curves in these figures are for the N(L)DPE model with 
/i^ = 0.3 and xjp/p^p < 0.1). Since D0 does not explicitly measure xjp/p^p, one explanation 
for their i^T-spectra is the F(IS)DPE model with the higher upper limits on xjp/p^p, say 
between 0.05 — 0.1. In such a case, recall that the interpretation of pomeron dominated 
exchange enters into question. For D018OO larger upper limits on xjp/p^p are inconsistent 
with the largest DPE dijet Et that is found, ~ 52 GeV. This maximum is consistent with 



an upper bound of only xjp/p^p < 0.057. On the other hand, for D063O, note from |]45|j48[ 
that dijets are produced up to Et ~ 25 GeV. If background and resolution effects or any 
other experimental complication can be ruled out, this suggests for D063O xjp/p^p ~ 0.08. 

For the y_|_-spectra, the CDF results unambiguously contradict interpretation as 
N(L)DPE dominated. The CDF y^-spectra is very broad, ranging as —2 ~ ?/+ ~ 1.5. 
From Fig. |^a, this is similar to the F(IS)DPE cases and is completely at odds with the 
N(L)DPE cases. Recall for the CDF cuts, the narrower ?/_|_-spectra for the N(L)DPE model 
is intrinsic to its lossless nature. Thus the experimental |/+-spectra is strong indication that 
the F(IS)DPE process dominates in the CDF case. This interpretation also is consistent for 
the ii^T-spectra in the most recent CDF report [53]. However, for the -Ey-spectra from the 



earlier CDF reports [^,^,^,0], ACTW SG is the most consistent model. Finally, one can 
not exclude the possibility that the N(L)DPE process gives a nondominant but measurable 
effect. This would help flatten any of the F(IS)DPE models in Fig. ^. However an admix- 
ture of the N(L)DPE process also will imply greater enhancement of the ?/+-spectra in the 
region — 1 ~ ?/+ ~ 0. 

Finally, total cross sections can be compared. CDF has preliminarily reported a DPE 
dijet cross section of (JTfPE ~ 13.6 nb |^. D0 has given no preliminary cross section. 



although for D018OO they have estimated ct^pe ~ 10 pb [0. In this case, the CDF total 



cross section is about 1000 times larger than D018OO. The model that most closely obtains 
this factor difference is F(IS)DPE ACTW D, which predicts the ratio to be still a factor 5 
smaller. As another case, for ACTW SG, it predicts the CDF total cross section should be 
only a factor ~ 50 larger than D018OO. Furthermore, the predicted ratio between the CDF 
and D018OO total cross sections decreases if the limits on xjp/p^p for D018OO are increased. 
Thus there appears to be some discrepancy between the experimental cross sections and 
those predicted by all the F (IS) DPE models. Also, the N(L)DPE models do very poorly in 
predicting the observed ratio. For the largest ratio predicted by these models, the CDF total 
cross section would be only a factor 7 larger than D018OO. However for the N(L)DPE case, 
modifications of the proton form factor and two-gluon pomeron model, which are discussed 
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in the Conclusion, may change these predictions by an order of magnitude but not more. 



V. CONCLUSION 

This paper has examined the factorized Ingelman-Schlein and nonfactorized lossless DPE 
processes of dijet production and computed predictions for the cuts of CDF, D0 and rep- 
resentative cuts of LHC. Two qualitative features emerge from our calculations which are 
model independent. They are reflections of the lossless kinematics of jet production in the 
N(L)DPE process, which requires that all the momentum carried by the pomerons will go 
into the hard event. The first of these qualitative differences between the N(L)DPE and 
F(IS)DPE processes emerges in the CDF ^/^-distributions in Figs. |^a and |]b. For this 
case, it is evident that for the N(L)DPE process, the distribution is considerably localized 
to within one unit of y^ rapidity whereas for the F(IS)DPE process, its distribution is con- 
siderably broader. The difference in the shapes of the distributions for these two processes 
is due to the intrinsic differences in the hard kinematics for the two processes and thus is 
expected to survive any nonperturbative modifications of the basic models. A second quali- 
tative difference emerges from the lower bound dependence of xjp for the N(L)DPE process. 
In our calculations, this difference was explicitly seen only amongst the set of LHC cuts, 
although the basic feature is general. Because of the lossless kinematics, if the lower bound 
on xjp is sufficiently large, low- ii^T jets will be prohibited from productions. An example of 
this is seen in the LHC ii^T-spectra, Figs. ||c and |]d, where for the 1' 2' cuts {xp > 0.00006) 
no truncation of low- ii^T jets is seen, whereas for the 1,2 cuts, where the lower bound on xjp 
is larger, xp > 0.002, a truncation of low-i^j^ jets occurs. 

Two additional qualitative features were found from our calculations which have some 
degree of model dependence, but for which we expect the general trends to sustain. The 
first of these, which appears for all the cuts, is that the Et spectra are fiatter (or "harder") 
for comparable N(L)DPE versus F(IS)DPE processes. This difference occurs due to the 
additional /^-dependence in the F(IS)DPE process that arises from the pomeron parton 
distribution functions. At fixed ^/s, larger Et implies larger (3, and the parton distribution 
functions fall, generally quickly, with increasing (3. There are modifications ||3^ to the basic 
N(L)DPE model, which will be discussed below, that can introduce a similar type of /?- 
dependence. As such, this difference between the N(L)DPE and F(IS)DPE processes has 
some model dependence. The second qualitative difference is in regards to the total cross 
sections, in particular the ratio of the CDF to D018OO total cross sections. The N(L)DPE 
process generally has a much smaller ratio for o'cdf/o'd0isoo ~ 4 compared to the F(IS)DPE 
process for which this ratio is typically ~ 100. However, any modeling that relies on the 



hard kinematics can alter this result, such as the /^-dependent modeling in [^ . 

In regards to our comparisons with the preliminary experimental data, we feel no final 
conclusions are possible, but the trends in the data show slight preference for domination 
by the F(IS)DPE process. From comparison of our models with either the recent ||5^ 



or earher [^,^,^,|^ CDF preliminary data, the mean rapidity (y+) spectra provide the 
most suggestive evidence that this data is dominated by the F(IS)DPE process. On the 
other hand, the D0 data shows no greater preference for either of the two mechanisms. 
Furthermore, none of the theoretical models examined here are completely consistent with 
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all the data. In particular, attempting to conclude that the data is dominated by the 
F(IS)DPE process is inconsistent with the harder i^^-spectra found in the D0 data and the 
earlier CDF data, since such spectra do not readily agree with the best fit ACTW F(IS)DPE 
model. Also, and perhaps most interesting, the experimental ratio of the CDF to D018OO 
total cross sections is a factor 5 — 10 larger than predicted by any of the models examined 
here. In this respect, it is worth noting that from Table 1, ratios between the D063O and 
D018OO cross sections also can be obtained, and once D0 complete their analysis of their 
preliminary data, this may be the next test between the predicted ratios and experiment. 
Finally, although it appears there are more indications that the data is dominated by the 
F(IS)DPE process, for the N(L)DPE process, neither data nor our analysis is sufficiently 
precise to rule out a subdominant component. All these issues are important to resolve as 
further experimental data becomes available. For the time being, the limitations of both 
the CDF and D0 experiments for this process exclude any final conclusions to be drawn. 

Our study emphasized the hard physics, but made minimal attempt at modeling most of 
the nonperturbative soft physics. For example, general belief is hard, pure hadron induced, 
diffractive processes are subject to a weakly i/i-dependent suppression factor [ p2| , pl| , which 



represents the probability for the rapidity gap(s) not to be filled by extra exchanges of 
pomerons and gluons between the particles in the model which have very different rapidities. 



Such absorptive corrections potentially can be treated by the methods developed in |32]. In 



39| , they used the approach of ||3^ and [^ to estimate the absorptive correction factor for 



the DPE dijet process and at ^/s = 1800 GeV found it to be ~ 0.06. For our calculations, 
this factor is meant to multiply the results in Table 1. 

For the N(L)DPE model, also not treated here is modifications of the two gluon exchange 
model with LLA ladder evolution of gluons. A model for this was given in [^]. Their 
model amounts to including in the N(L)DPE dijet amphtude two factors of gluon densities 
evaluated at xjp/p and xjp/p. This modification has interesting consequences. For example, 
maintaining the same procedure to normalize the DPE dijet amplitude with respect to the 
elastic pp cross section, the N(L)DPE cross section with this modification then will be 1 — 2 
orders of magnitude smaller. This arises because the normalization constant is fixed for 
elastic scattering process, where x^^^p ~ l±/v^ is small since 1^, the transverse momentum 
of the outgoing hadrons, is ~ 1 GeV. For example at ^/s = 1800 GeV, x^y^^p ~ 0.0006. 
On the other hand, for the N(L)DPE process the typical xjp/p^p are ~ 0.01. These are 
much larger, which implies a decrease of the gluon densities relative to the elastic scattering 
situation. This provides an additional suppression compared to the same model with no 



gluon densities. In [^, the exclusive double diffractive dijet model is the same as our 
N(L)DPE model here, except for the inclusion of gluon densities in their model. It is due to 
the effect of these densities that their model predicts cross sections 1 — 2 orders of magnitude 
smaller than ours. 

Another effect of the gluon densities is the N(L)DPE i^-^-spectra will fall faster with 
Et- The reasons for this are the same as explained in Sect. |^ for the F(IS)DPE models. 
In short, higher Et requires larger xjp/p^p which in turn implies smaller gluon densities. 
This effect will be less pronounced for the N(L)DPE models versus the F(IS)DPE models, 
since the argument of the gluon densities in the former always remain small, ~ x^'^^ ~ 0.1 
whereas in the latter it ranges up to ~ 1, where the maximum diminution of the gluon 
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densities occurs. 

In summary, this paper examined two very different types of diffractive mechanisms 
for DPE dijet production, the factorized Ingelman-Schlein and nonfactorized lossless mech- 
anisms. In the spirit of Regge physics, both mechanisms were termed double pomeron 
exchange. Some of the differences between the two processes have been elucidated here, 
which will help in interpreting experiment. The F(IS)DPE model appears best to represent 
the present experimental data. However, inconsistencies still remain that need to be sorted- 
out by both theory and experiment before final conclusions can be made. Although the 
N(L)DPE process does not appear to dominate the cross section, a subdominant component 
of it can not be excluded with present information. The cleanliness of the final state in 
this process, two outgoing hadrons plus a hard event, suffices justification to search for it. 
For new particle search experiments, this process could permit the ultimate measurement, 
although its diminutive cross section precludes it from being the ideal measurement. Thus 
as should be the case with any good tale, the hard double pomeron exchange story is filled 
with uncertainties and confiicts. Furthermore, the hope vested in the N(L)DPE process 
someday may vindicate a familiar moral, that the best things don't come easy. 
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Table 1: Double Pomeron Exchange (DPE) dijet total cross sections for the nonfactorized 
(N) and factorized (F) models. 
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FIG. 1. Double pomeron exchange (DPE) to two jets. 




FIG. 2. Amplitude for inclusive two-jet production. 




FIG. 3. Factorized (Ingelman-Schlein) Double Pomeron Exchange (F(IS)DPE) amplitude with 
two jets produced. 








FIG. 4. Our model of the Nonfactorizing (Lossless) Double Pomeron Exchange (N(L)DPE) 
amplitude with two gluon jets produced. 
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FIG. 5. Double pomeron exchange (DPE) dijet production-CDF cuts, (a) Mean rapidity 
spectra y+ = (yi + y2)/2. (b) Ratio of mean rapidity spectra between the DPE and correspond- 
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FIG. 6. Double pomeron exchange (DPE) dijet production - D018OO cuts, (a) Mean rapidity 
spectra ^+ = (^i + ^2)/2. (b) Ratio of mean rapidity spectra between the DPE and correspond- 
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FIG. 7. Double pomeron exchange (DPE) dijet production - D063O cuts, (a) Mean rapidity 
spectra y+ = (%j\ + y2)/2. (b) Ratio of mean rapidity spectra between the DPE and correspond- 
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FIG. 8. Double pomeron exchange (DPE) dijet production - LHC cuts, (a) Mean rapidity 
spectra y+ = (yi + y2)/2. (b) Ratio of mean rapidity spectra between the DPE and correspond- 



